Fighting Multidrug Resistance with Ruthenium-Cyclopentadienyl Compounds: Unveiling the Mechanism of P-gp Inhibition.

The search for more effective and selective drugs to overcome cancer multidrug resistance is urgent. As such, a new series of ruthenium-cyclopentadienyl ("RuCp") compounds with the general formula [Ru(η5-C5H4R)(4,4'-R'-2,2'-bipy)(PPh3)] were prepared and fully characterized. All compounds were evaluated toward non-small cell lung cancer cells with different degrees of cisplatin sensitivity (A549, NCI-H2228, Calu-3, and NCI-H1975), showing better cytotoxicity than the first-line chemotherapeutic drug cisplatin. Compounds 2 and 3 (R' = -OCH3; R = CHO (2) or CH2OH (3)) further inhibited the activity of P-gp and MRP1 efflux pumps by impairing their catalytic activity. Molecular docking calculations identified the R-site P-gp pocket as the preferred one, which was further validated using site-directed mutagenesis experiments in P-gp. Altogether, our results unveil the first direct evidence of the interaction between P-gp and "RuCp" compounds in the modulation of P-gp activity and establish them as valuable candidates to circumvent cancer MDR.

Histidine Protonation and Conformational Switching in Diphtheria Toxin Translocation Domain.

Protonation of key histidine residues has been long implicated in the acid-mediated cellular action of the diphtheria toxin translocation (T-) domain, responsible for the delivery of the catalytic domain into the cell. Here, we use a combination of computational (constant-pH Molecular Dynamics simulations) and experimental (NMR, circular dichroism, and fluorescence spectroscopy along with the X-ray crystallography) approaches to characterize the initial stages of conformational change happening in solution in the wild-type T-domain and in the H223Q/H257Q double mutant. This replacement suppresses the acid-induced transition, resulting in the retention of a more stable protein structure in solutions at pH 5.5 and, consequently, in reduced membrane-disrupting activity. Here, for the first time, we report the pKa values of the histidine residues of the T-domain, measured by NMR-monitored pH titrations. Most peaks in the histidine side chain spectral region are titrated with pKas ranging from 6.2 to 6.8. However, the two most up-field peaks display little change down to pH 6, which is a limiting pH for this protein in solution at concentrations required for NMR. These peaks are absent in the double mutant, suggesting they belong to H223 and H257. The constant-pH simulations indicate that for the T-domain in solution, the pKa values for histidine residues range from 3.0 to 6.5, with those most difficult to protonate being H251 and H257. Taken together, our experimental and computational data demonstrate that previously suggested cooperative protonation of all six histidines in the T-domain does not occur.

Interfacial Dynamics and Growth Modes of ß2-Microglobulin Dimers.

Protein aggregation is a complex process, strongly dependent on environmental conditions and highly structurally heterogeneous, both at the final level of fibril structure and intermediate level of oligomerization. Since the first step in aggregation is the formation of a dimer, it is important to clarify how certain properties of the latter (e.g., stability or interface geometry) may play a role in self-association. Here, we report a simple model that represents the dimer's interfacial region by two angles and combine it with a simple computational method to investigate how modulations of the interfacial region occurring on the ns-µs time scale change the dimer's growth mode. To illustrate the proposed methodology, we consider 15 different dimer configurations of the ß2m D76N mutant protein equilibrated with long Molecular Dynamics simulations and identify which interfaces lead to limited and unlimited growth modes, having, therefore, different aggregation profiles. We found that despite the highly dynamic nature of the starting configurations, most polymeric growth modes tend to be conserved within the studied time scale. The proposed methodology performs remarkably well taking into consideration the nonspherical morphology of the ß2m dimers, which exhibit unstructured termini detached from the protein's core, and the relatively weak binding affinities of their interfaces, which are stabilized by nonspecific apolar interactions. The proposed methodology is general and can be applied to any protein for which a dimer structure has been experimentally determined or computationally predicted.

Novel US-CpHMD Protocol to Study the Protonation-Dependent Mechanism of the ATP/ADP Carrier.

We have designed a protocol combining constant-pH molecular dynamics (CpHMD) simulations with an umbrella sampling (US) scheme (US-CpHMD) to study the mechanism of ADP/ATP transport (import and export) by their inner mitochondrial membrane carrier protein [ADP/ATP carrier (AAC)]. The US scheme helped overcome the limitations of sampling the slow kinetics involved in these substrates' transport, while CpHMD simulations provided an unprecedented realism by correctly capturing the associated protonation changes. The import of anionic substrates along the mitochondrial membrane has a strong energetic disadvantage due to a smaller substrate concentration and an unfavorable membrane potential. These limitations may have created an evolutionary pressure on AAC to develop specific features benefiting the import of ADP. In our work, the potential of mean force profiles showed a clear selectivity in the import of ADP compared to ATP, while in the export, no selectivity was observed. We also observed that AAC sequestered both substrates at longer distances in the import compared to the export process. Furthermore, only in the import process do we observe transient protonation of both substrates when going through the AAC cavity, which is an important advantage to counteract the unfavorable mitochondrial membrane potential. Finally, we observed a substrate-induced disruption of the matrix salt-bridge network, which can promote the conformational transition (from the C- to M-state) required to complete the import process. This work unraveled several important structural features where the complex electrostatic interactions were pivotal to interpreting the protein function and illustrated the potential of applying the US-CpHMD protocol to other transport processes involving membrane proteins.

Predicting stable binding modes from simulated dimers of the D76N mutant of ß 2-microglobulin.

The D76N mutant of the ß 2 m protein is a biologically motivated model system to study protein aggregation. There is strong experimental evidence, supported by molecular simulations, that D76N populates a highly dynamic conformation (which we originally named I 2 ) that exposes aggregation-prone patches as a result of the detachment of the two terminal regions. Here, we use Molecular Dynamics simulations to study the stability of an ensemble of dimers of I 2 generated via protein-protein docking. MM-PBSA calculations indicate that within the ensemble of investigated dimers the major contribution to interface stabilization at physiological pH comes from hydrophobic interactions between apolar residues. Our structural analysis also reveals that the interfacial region associated with the most stable binding modes are particularly rich in residues pertaining to both the N- and C-terminus, as well residues from the BC- and DE-loops. On the other hand, the less stable interfaces are stabilized by intermolecular interactions involving residues from the CD- and EF-loops. By focusing on the most stable binding modes, we used a simple geometric rule to propagate the corresponding dimer interfaces. We found that, in the absence of any kind of structural rearrangement occurring at an early stage of the oligomerization pathway, some interfaces drive a self-limited growth process, while others can be propagated indefinitely allowing the formation of long, polymerized chains. In particular, the interfacial region of the most stable binding mode reported here falls in the class of self-limited growth.

pK a Calculations in Membrane Proteins from Molecular Dynamics Simulations.

The conformational changes of membrane proteins are crucial to their function and usually lead to fluctuations in the electrostatic environment of the protein surface. A very effective way to quantify these changes is by calculating the pK a values of the protein's titratable residues, which can be regarded as electrostatic probes. To achieve this, we need to take advantage of the fast and reliable pK a calculators developed for globular proteins and adapt them to include the explicit effects of membranes. Here, we provide a detailed linear response approximation protocol that uses our own software (PypKa) to calculate reliable pK a values from short MD simulations of membrane proteins.

Improved GROMOS 54A7 Charge Sets for Phosphorylated Tyr, Ser, and Thr to Deal with pH-Dependent Binding Phenomena.

Phosphorylation is a ubiquitous post-translational modification in proteins, and the phosphate group is present constitutively or transiently in most biological building blocks. These phosphorylated biomolecules are involved in many high-affinity binding/unbinding events that rely predominantly on electrostatic interactions. To build accurate models of these molecules, we need an improved description of the atomic partial charges for all relevant protonation states. In this work, we showed that the commonly used protocols to derive atomic partial charges using well-solvated molecules are inadequate to model the protonation equilibria in binding events. We introduced a protocol based on PB/MC calculations with a single representative conformation (of both protonation states) and used the resulting pKa estimations to help manually curate the atomic partial charges. The final charge set, which is fully compatible with the GROMOS 54A7 force field, proved to be very effective in modeling the protonation equilibrium in different phosphorylated peptides in the free (tetrapeptides, pentapeptides, and pY1021) and protein-complexed forms (pY1021/PLC-γ1 complex). This was particularly important in the case of the pY1021 bound to the SH2 domain of PLC-γ1, where only our curated charge set captured the correct protonation equilibrium at the neutral to slightly acidic pH range. The binding/unbinding phenomena in that pH range are biologically relevant, and to improve our models, we need to go beyond the commonly used protocols and obtain revised force field parameters for these molecules.